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Abstract 

We apply the concepts of utility based pricing and hedging of derivatives in 
stochastic volatility markets and introduce a new class of "reciprocal affine" models 
for which the indifference price and optimal hedge portfolio for pure volatility claims 
are efficiently computable. We obtain a general formula for the market price of 
volatility risk in these models and calculate it explicitly for the case of an exponential 
utility. 
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1 Introduction 



The purpose of this paper is to provide a concrete application of utility based pricing and 
hedging for markets with untraded risk factors. For incomplete markets of this type, they 
provide an economically justified framework for pricing derivatives where the traditional 
paradigms of replication and risk neutral valuation fail to produce any identifiable price. 
For instance, this framework is well suited for pricing in insurance markets (see Young 
and Zariphopoulou 2002) as well as market models with stochastic volatility. 

The method is based on solutions of the so called optimal hedging problem (see Cvi- 
tanic, Schachermayer and Wang 2001; Owen 2002), which is a generalization of Merton's 
optimal investment problem, whose solution has a distinguished history including key 
works such as Karatzas, Lehoczky, Shreve and Xu (1991) , and Kramkov and Schacher- 
mayer (1999) and Schachermayer (2001). Characteristically, these papers present abstract 
existence and uniqueness results for the optimal portfolios, but do not address the ques- 
tion of actually constructing or systematically approximating them given specific market 
models. However, because of the variety and richness offered by general incomplete mar- 
kets, a full understanding of the appropriateness of this pricing concept can only achieved 
by applying it to an increasing number of concrete examples. 

The valuation of derivatives in markets with stochastic volatility typically depends 
on the choice of an appropriate dynamics for the volatility process a t , the specification 
of a functional form for the volatility risk premium and the use of asymptotic filtering 
techniques to estimate the model parameters and the current level of the unobservable 
process a t in terms of the time series of observed stock prices S t (see Fouque, Papanicolaou 
and Sircar 2000; Griinbichler and Longstaff 1996; Heston 1993; Stein and Stein 1991). 
While we do not address estimation and calibration issues, we show how the utility based 
framework offers insight for the choice of stochastic volatility models and provides a 
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specification of volatility risk premia. 

We consider contingent claims on two factor stochastic volatility models of the form 



dS t = SMt,Y t )dt + a{t,Y t )dW}], 

dY t = a(t, Y t )dt + b(t, Y t ) [pdWl + y/l - p 2 dW?}, (1) 

with initial values So, Y > 0, deterministic functions fi, a, a, b satisfying the regularity 
and growth conditions necessary to ensure the existence and uniqueness of a solution to 
(|T|). independent one dimensional P-Brownian motions and W 2 , and finally constant 
correlation \p\ < 1. 

In section 2, we review the concepts of certainty equivalent and indifference price of a 
claim as introduced by Hodges and Neuberger (1989) (according to Becherer 2001) ) for 
a general utility function and then specialize to an exponential utility U (x) = —e~ ix with 
risk aversion parameter 7 > 0. We then show how the utility based approach uniquely 
specifies the volatility risk premium associated with each contingent claim B(St,Yt), as 
well as a claim independent market price for volatility risk which is related to the solution 
of Merton's problem for a stochastic volatility market. 

We use the theory of dynamic programming in section 3 to derive the nonlinear PDE 
satisfied by the certainty equivalent of a general contingent claim on the market specified 
by (0). When the contingent claim is a pure volatility claim (such as in Brenner and Galai 
1989; Whaley 1993), we will see in section 4 how a change of variables leads to a linear 
second order PDE, for which a Feynman-Kac representation of the solution is available. 

In section 5, we explore the striking resemblance of this formula to the framework of 
interest rate derivatives. We introduce the reciprocal affine stochastic volatility model, 
where the "spot rate" R t oc a~j~ 2 is taken to be a CIR process. 

In section 6 we use a transform analysis technique developed by Duffle, Pan and 
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Singleton (2000) to obtain tractable expressions for the indifference price and the holdings 
in both the hedging and Merton portfolios for stochastic volatility. The choice of a 
reciprocal affine model for the squared volatility process fully reveals its intrinsic elegance 
in this section when we compute its market price for volatility risk explicitly. In section 
7 we benefit from the efficiently implementable formulas derived previously and present 
numerical results for volatility put options and the difference of volatility call options, 
discussing several qualitative features of the indifference prices and the hedging strategies 
for such instruments. In particular, for 7 not too large, we observe clearly that pricing 
and hedging are approximately independent of the risk aversion. 

2 Utility based pricing 
2.1 General setup 

The central definition is that of an optimal hedging portfolio. It corresponds to the strategy 
followed by an investor with initial wealth x who, when faced with a (discounted) financial 
liability B maturing at a future time T, tries to solve the stochastic control problem 

u(x) = sup E [U (X T - B) \X = x] , (2) 

A 

where X T is the discounted terminal wealth obtained when investing H t S t dollars on the 
risky asset and r) t C t dollars in a riskless cash account with value C t initialized at Co = 1 
and governed by 

dC t = r t C t dt. (3) 

The utility function U : R — * R U {—00} is assumed to be a concave, strictly increasing 
and differentiable function and the domain of optimization A is restricted to self-financing 
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portfolios, that is, to wealth processes satisfying 

C t X t :=H t S t + Vt C t = x+ [ H u dS u + [ Vu dC u . (4) 

Jo Jo 

In addition to the self-financing condition, economic reasoning imposes further restrictions 
on the class A of admissible portfolios. For instance, one might require that the wealth 
process X t be uniformly bounded from below, in order to rule out "doubling strategies" 
(following Harrison and Pliska 1981). Finally, the liability B is assumed to be a random 
variable of the form B = B(St, Yt), for some Borel function B : IR+ — > R. 

It follows from Ito's formula and the definitions (0) and Q that the market model in 
terms of the discounted price process S t = S t /C t is 

dS t = St[(ti(t,Yt)-r)dt + <r(t,Y t )dW}], 

dY t = a(t, Y t )dt + b(t, Y t ) [pdW t l + p 2 dW?\. (5) 

We have taken r t = r to be constant for simplicity; the case of either a deterministic 
rate r t = r(t) or a stochastic rate of the form r t = r(t, St, Y t ) are immediate but artificial 
generalizations. The really relevant case of stochastic rates driven by a different stochastic 
factor, possibly correlated with St and Y t , is not pertinent to the scope of this paper. 

Using the self- financing condition (|4"J) . we immediately obtain that the discounted 
wealth process satisfies 

dX t = H t dSt = H t S t [(jji(t, Y t ) - r)dt + a(t, Y t )dW}], (6) 

so that the only relevant control in (J2J is H t , with the holdings in the cash account being 
determined by rjt = X t — H t S t . 
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For Markovian markets such as (jSj) and claims of the form B = B(St,Yt), we can 
embed the optimal hedging problem (j2J) into the larger class of optimization problems 
defined by 

u(t,x,s,y)= sup E t ^ y [U(X T -B(S T ,Y T ))\X t = x], (7) 
HeAt 

for t G (0, T), where x G 1R denotes some arbitrary level of wealth, At denotes admissible 
portfolios starting at time t and E ttS , y [-\ denotes expectation with respect to the joint 
probability law at time t of the processes S u , Y u satisfying (jSJ) for u > t, with initial 
condition S t = s and Y t = y. 

2.2 The certainty equivalent and the indifference price 

Suppose that (JJJ) has an optimizer Hf , that is, assume that 

u(t, x, s, y) = E tiS>y [U(x + (H B ■ S)f - B(S T , Y T ))], 
where we use the notation 

(H-S)J := £ H u dS u . 

Define the certainty equivalent for the claim B at time t as the process cf = c B (t, x, s, y) 
satisfying the equation 

U(x - cf )) = E t>s , y [U(x + (H B ■ S)f - B(S T , Y T ))]. (8) 

In other words, it corresponds to the amount which when subtracted from the wealth x 
at time t produces the same deterministic value for the utility as the optimal expected 
utility of terminal wealth starting at x and facing the claim B at terminal time T. If we 
set 5 = 0, then the optimal hedging problem becomes the Merton optimal investment 



6 



problem and we denote the certainty equivalent by c° = c°(t, x, s, y). 

An agent with utility U and wealth x at time t G (0, T) will charge a premium for 
issuing a liability B maturing at T . The indifference price for the claim B is defined to 
be the premium that makes the agent indifferent between making the deal or not, that 
is, the unique solution n B = n B (t,x, s,y) (if it exists) to the equation 

sup E t , s , y [U(x + (H ■ S)J] = sup E t>a>y [U(x + n B + (H ■ S)J - B(S T , Y T )]. (9) 
He At HeAt 

From the definition of the certainty equivalent, we see that this equation is equivalent to 

U(x-c° t ) = U(x + 7r B -c?), (10) 



so that the indifference price is given by 

7r s (t, x, s, y) = c B (t, x + ir B (t, x, s, y), s, y) - c°(t, x, s, y). (11) 



The Davis price (see Davis 1997) of the contingent claim B is defined as the linearization 
of the indifference price, that is, 



b . = d(^ B ) 

^ Davis • 



Javis ■ j 



(12) 



As we see, establishing the existence of an indifference price is tantamount to solving 
both the optimal hedging problem for the claim B and the Merton investment problem. 
As such, it depends on the choice of the utility function and the appropriate notion of 
admissible portfolios, as well as integrability conditions on the claim B itself. 
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For the rest of this paper we concentrate on exponential utilities of the form 



U(x) = -e" 7X , (13) 

where 7 > is the risk-aversion parameter. The advantage of using an exponential utility 
is that we can factorize the value function u(t, x, s, y) in (J7J) as 



u(t,x,s,y) = sup E tjS! y 
He At 



-e 1X inf E t ,s, y 



e -y(x+(H^S)f-B(S T ,Y T )) 

e -~,((H-S)T-B(S T ,Yt)) 



=: U(x)v{t,s,y). (14) 

To define the domain of optimization, let M a (S) and M e (S) denote respectively the 
sets of absolutely continuous and equivalent (local) martingale measures for S and let 
M-^(S) denote the set of measures Q G A4 a (S) with finite relative entropy with respect 
to P. For absence of arbitrage, we assume that M e H M f ^0. Now let L(S) denote the 
set of predictable S-integrable processes. For definiteness, we take the set of admissible 
portfolios to be 

A={H e L(S) : (H ■ S) t is a Q-martingale for all Q G Mf), (15) 

but note that alternative choices for A also yield the existence and uniqueness of optimal 
hedging portfolio for the exponential utility. It follows from the abstract convex duality 
solution of the optimal hedging problem that, for claims B satisfying certain integrabil- 
ity conditions, the certainty equivalent derived using the exponential utility exists as a 
well defined process (see Becherer 2001, 2004; Delbaen, Grandits, Rheinlander, Samperi, 
Schweizer and Strieker 2002; Kabanov and Strieker 2002; Owen 2002). 
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It follows directly from (jHJ) that the certainty equivalent is wealth independent and 
given by 

(16) 



c (t,s,y) = - logv(t,s,y). 

7 

Analogously, setting B = gives the certainty equivalent for the Merton problem with 
exponential utility as 

(17) 



c°(t,s,y) = - logv°(t,s,y), 

7 



where v°(t, s, y) — inf E t> 



HeAt 



s.y 



Thus the indifference price process for the claim 



B obtained from an exponential utility is given by 



n B (t, s, y) = c B (t, 8 , y) - c°(t, s,y) = - log VU ' " ^ 



7 v°(t,s,y)' 



2.3 Martingale Pricing 

The solution of the hedging problem (J2J) through convex duality involves finding a martin- 
gale measure Q B minimizing the expected value of the Legendre transform of the utility 
function U over the set Ai a (S) of absolutely continuous martingale measures for S. The 
optimal wealth X B = x + (H B ■ S)t and the optimal martingale measure Q B are then 
related by the fundamental equation 



U'{Xf -B)=( 



dQ B 
dP 



(19) 



where £ = u'(x) (see Owen 2002). Letting 



A 



B 



dQ B 
dP 



-E t [U'(X B -B)) 



(20) 



9 



be the density process for the measure Q B , we define the utility based price of risk asso- 
ciated with the claim B as the vector process \f = ((A B )*, (\ B )t) satisfying 

| = -[(« + («. (21) 

For the case of exponential utility this reduces to 



Af = -l Et [U(X t B + (H B -S)J-B] 
= -^U(X B -cf) 



2 - 7 (Xf-cf) 
e - 7 (x- C S) ' 



(22) 



where we have used the definition of the certainty equivalent process cf and the fact that 
£ = u'{x) = U'(x — cf ). Applying Ito's lemma to the martingale Af above, we easily 
obtain that 



(\ B )l = ^-d^aS-bpdyC 3 }, (23) 
{\ B f t = - 7 ^c B yr^7. (24) 

Both the optimal martingale measure Q B and the utility based price of risk Af are 
specifically related to the claim B and therefore do not constitute a direct generalization 
of the paradigm of pricing by expectation with respect to a risk adjusted measure valid 
for all claims. For instance, as we will later observe in our numerical examples, this is 
reflected in the fact that the indifference price n B is not linear in the claim B, with the 
obvious effect that Q B fails to be a pricing measure even for multiples of B, let alone for 
other unrelated claims. For the Davis price, however, matters are slightly more familiar, 
as the following argument shows. 
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By differentiating the identity 



U(x - cf) = E t [U{x + (H £B ■ S)J - eB)} 



at e = we obtain Davis's formula 



dc £ t 



Davis 



E t [U'(x + (H°-S)f)B] 
U'{x-c° t ) 

E?[B), 



(25) 



where we have used (|19j) for Merton's problem. We see that the Davis price can be 
calculated as an expectation with respect to the risk adjusted measure Q obtained as 
the dual solution for the optimal investment problem according to the utility function U. 
This measure induces a utility based market price of risk valid for all claims, which in the 
case of an exponential utility is obtained from by setting B = 0. 

3 The PDE for the certainty equivalent 

Consider the continuous Markovian market specified by (0), an exponential utility function 
of the form (JT3j) and a claim of the form B = B(St, Yt) satisfying appropriate integrability 
conditions. By the dynamic programming principle, the value function u(t, x, s, y) satisfies 
the Hamilton- Jacobi-Bellman equation 

1 1 

d t u + -s 2 a 2 d 2 ss u + bpsad 2 sy u + -b 2 d 2 yy u + s(ii-r)d s u + ad y u (26) 



+ max < -h 2 s 2 a 2 d 2 x u + bphsad 2 u + hs 2 a 2 d 2 s u + hs(/j, — r)d x u 




0. 
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with boundary condition u(T, z) = — e 1<yX B ( s > y " and the optimal portfolio process is 
given by 

Hf = h B {t,X t ,S t ,Y t ), (27) 

where h B is the optimizer of the expression above. Direct substitution of (|14|) into this 
HJB problem leads to an optimizer of the form H B = h B (t, s, y), with 



h B (t,s,y) = -— + 



1 d s v bp d y v (// — r) 



+ 



7 v 7scr v 



, ysa z 



(2f 



The partial differential equation satisfied by the optimal function v(t,s,y) is then 



d t v + - (s 2 a 2 d 2 ss v + 2bpsad 2 ys v + b 2 d 2 yy v) + 



bp(p - r) 



dyV 



-{bpdyV + sad s v) 2 + 



(/i-r) s 



0. 



(29) 



subject to the terminal condition v(T, s,y) = e 13 ^^. 

From (JTHj), we find that the certainty equivalent process c B (t, s,y) is a solution to the 
partial differential equation 



d t c B + ^( s V9 s 2 s c B + 2 S a6p^ C B + 6 2 ^c B ) + 



yy 



bp(p — r) 



a 



dyC 



%^ + ^(l-p 2 )(^) 2 



(30) 



with terminal condition c B (T,s,y) = B(s,y). The nonlinearity (d y v) 2 /v for the value 
function has been transformed into the quadratic nonlinearity (d y c B ) 2 for the certainty 
equivalent, while the nonlinearity associated with (d s v) 2 /v is absent. The partial differ- 
ential equation satisfied by c°(t,s,y), the certainty equivalent for Merton's problem, is 
identical to ([3*0)1. but with the terminal condition c°(T, s,y) = 0. 
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Using (j28|). the optimal portfolio can be obtained in terms of the certainty equivalent 
process by 

h B (t, s, y) = d s c B + ^d y c B + (31) 

From (|23|l . (|24jl we obtain that the exponential utility based price of risk associated with 
the claim B is 

mbu _ W,Yt)~r) , , 

0\t, Y t ) 

= -TH^v^?, (33) 

with the desirable property that the dependence on the claim B occurs only in its second 
component, while the component related to the first Brownian motion is formally identical 
to the market price of risk of a complete market. Setting B = in the equations above 
gives the market price of volatility risk induced by the exponential utility through the 
solution of Merton's optimal investment problem. 



4 Volatility claims 

The equations from the previous section simplify considerably when the claim is inde- 
pendent of the process St- Because of the exponential nature of the SDE for St in (JHJ), 
we obtain that for pure volatility claims of the form B = B(Yr), the equation for the 
certainty equivalent cf = c B (t,y) is reduced to 



d t c B + 



bp(n - r) 
a 



a 



dy^ + )fV - (J 0- + - P 2 )^) 2 = °' (34) 
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subject to the terminal condition c B (T,y) = B(y), whereas the optimal hedging portfolio 
is given by 



h B (t,s,y) = - 



^d y c B + 



a 



7<j 2 



Following Zariphopoulou (2001) we now use the transformation 



(35) 



c B (t,y) 



1 



7(1 -P 2 ) 



log /(£,?/), 



to reduce (|34|) to the linear parabolic final value problem 



bp(p - r) 



a 



2a 2 



(36) 



(37) 



f{T,y) 



Under the appropriate growth and boundedness assumptions on the coefficient func- 
tions fi, a, a and b, we can use the Feynman-Kac formula to represent the solution to the 
problem above as 



where we define 



f{t,y) = E t>v 



R(t,y) 



e - jf R(.s,Ys)ds e j(l-p 2 )B(Y T ) 



{l-p 2 ){p{t,y)-rf 



(38) 



(39) 



2a(t,yy 

and E tt y[-} denotes the expectation with respect to the probability law at time s = t of 
the solution to 



dY s 



b{jJL — r)p 



a 



ds + b 



pdW} + v 7 ! - P 2 dW t 



(40) 



for a pair of independent one dimensional P-Brownian motions Wf, Wj~, for a probability 
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measure P on (fl, J 7 , (J-t)te[o,Tl)- If we further required S to be a P martingale, the 
comparison with (|SJ) leads to the identification 



i 



dWl + Ajd* 



dWl 



(41) 



where 



A 



»(t, Y t ) - r 
°(t,Y t ) 



(42) 



The martingale measure P is a computational device in order to obtain the function 
f(t,y) through (|3"8)l. and should not be confused with either the optimal martingale 
measure Q B associated with the claim B or the risk adjusted measure Q giving rise to 
the Davis price as explained in section 2. 

5 Reciprocal affine models 

The first factor in the integrand of the Feynman-Kac representation (|38|) is reminiscent 
of a stochastic discount factor with a "spot rate" Rt implicitly related to our original 
non-traded factor Y t through Following this analogy, we see that when B = 0, (|S5)l 
is the formal equivalent of the price of a zero coupon bond. If in addition equation (jHHJ) 
can be inverted, then the second factor in ()38|) becomes a claim on the spot rate Rt and 
we can use the technology developed for interest rate derivatives in order to price it. 

To explore this analogy with interest rates further, for the remainder of the paper we 
take /i and r to be constants and a(t, Y t ) = s/Y t , so that ()39|) becomes 



Rt = R(t,Y t ) 



{l-p 2 )^-r? 
2Y t 



(43) 
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Affine models form a well-studied class of interest rate models (see Duffie, Pan and 
Singleton 2000), often leading to analytic expressions for quantities such as bond prices. 
We can carry the results from these models to our problem by hypothesizing that Rt 
follows an affine process and thus Y t is a reciprocal affine process. We illustrate the idea 
in the specific case of the Cox, Ingersoll and Ross (1985) model. Since our calculations 
are going to take place under the measure P, we specify the dynamics for R t as 



for constants a,K,[3 > with 4an > (3 2 . It follows from (j41j) that the dynamics of R t 
under the economic measure P is 



dR t = a(H - Rt)dt + (3\[R t pdW} + yj\ - p 2 dW 2 , 



(44) 



dRt = a(K- Rt)dt + /3^R t pdW] + y/l - p 2 dW? , 



(45) 



where a 



and an = an. We then obtain from the Ito formula that 



dY t 




Y t 2 + aY t dt 



pdWl + v 7 ! - P 2 dW? , 



(l-p 2 )(p-r) 2 
2itf 



dR\ 



(46) 



which by comparison with (JHJ) leads to 



ait, Y t ) 



aY t + 



2((3 2 - an) 




(47) 



(1 - p 2 )(p-r) 




(48) 
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The reciprocal of a CIR process, being positive and mean reverting, is a reasonable 
model for stochastic volatility. Yet the expressions above show that our formulation differs 
from popular choices for stochastic volatility, such as the Heston (1993) model, for which 
Y t itself follows a CIR process, or the Stein and Stein (1991) model, where a t follows an 
arithmetic Ornstein-Uhlenbeck process. 

6 Pricing and hedging formulas 

The indifference price and Davis price of a volatility claim under the reciprocal affine 
models of the previous section both require computation of expressions of the form 

(49) 

for functions g : IR + — > KL Provided its Fourier transform is well defined and invertible, 
we can express g as 

1 f°° 

9{R) = ^J e-* uR g(u)du, (50) 

where 

/oo 
e iuR g(R)dR. (51) 
■oo 

Using Fubini's theorem to interchange the order of integration, we have 

1 f°° 

I = I(R t , t,T) = —J *(u)g(u)du, (52) 

where \1/ can be computed following Duffle, Pan and Singleton (2000) as 

*(u)=*(u,Rt,t,T) := E t 

= exp[M(u,t,T) + N(u,t,T)R t \. (53) 



I:=E t 



- J 4 T Rsds 



g(Ri 



■ J ( T R s ds —iuRt 
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Here 



N(u) = N(u,t,T) 
M(u) = M(u,t,T) 



(b 2 + iu)bx - (61 + w)6 2 e A( *" T) 
(6 2 + iuj - (61 + iu)e A (*- T ) ' 
— 2aK , fb 2 + iu\ , . mN 



(54) 



with 6 2 > 61 being the two roots of a; 2 — — and A = a/« 2 + 2/3 2 . 

Setting <7(i2r) = e -r^-^)B(Br) ; we obtain from (HHJ), (J2S1) and (03) that the indifference 
price of the volatility claim B = B(Rt) is simply 



7T 



B 



-log 

7 



7(1 -P 2 ) 



log 



I(Rt,t,T) 



(55) 



From (|35|1. ()36|) and ()38|h the number of shares of stock to be held in order to optimally 
hedge against the claim B is 



h B (t,y) 



1 

7s 



6p d log/ (/i — r) 



7(1-P 2 )v^ <%/ 



71/ 



J_(/x--r) 
7-5 2/ 



/5p n,*(«)iV(«^(tt)dti 



^(1-p 2 ) r M *(«)s(«)<i 



+ 1 



whereas the number of shares held in the Merton portfolio is 



(56) 



h°(t,y) 



bp <91og*(0) (fi-r) 



s [i(l-p 2 )y/y dy 



1 (l*-r) 
is y 



(3p 



\/2(l-p 2 ) 



N(0) + 1 



(57) 



Finally, using (J32j) with 5 = 0, the reader is invited to perform an amusing calculation 
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leading to the conclusion that the market price of volatility risk induced by an exponential 
utility in the reciprocal affine model is 



V Xt ) 



( 



\ 



(p-r) 



(58) 



This remarkable functional form for the market price of volatility risk has several pleasant 
implications. We first observe that in the regime of a sufficiently remote time to maturity 
both of its components are constant multiples of the usual market price of risk obtained for 
complete markets with a time varying volatility a t = y/Y t . Secondly, for typical parameter 
values, the first component is about two orders of magnitude bigger than the second one, 
implying that most of the volatility risk is already accounted by the dependence of the 
model on the first Brownian motion W}. Finally, it is now easy to verify that our process Y t 
is the reciprocal of a CIR process under both the economic measure P and the exponential 
risk adjusted measure Q, a consistency result which is far from obviously true for other 
stochastic volatility models. 



7 Numerical results 

Numerical computation of prices and hedge amounts is now highly efficient using the fast 
Fourier transform. We found that with a discrete Fourier lattice of 2 12 points, we could 
compute hundreds of prices and hedging amounts per second with accuracy better than 
1% on a desktop PC. This speed enabled us to easily survey dependences of prices and 
hedging in all parameter values. 

We illustrate the range of possibilities for model parameters fixed at reasonable values: 
p = 0.5, a = 5,/3 = 0.04, k = 0.001, p = 0.04, r = 0.02 and initial squared volatility 
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ranging in the interval [0,0.5]. With these parameters the squared volatility process has 
a mean reversion time of approximately two months and an equilibrium distribution with 
expected value approximately 40%. 

Figures 1, 2 and 3 present plots of the indifference price of a volatility put option with 
payoff (K — (Jt) + with strike K = 0.15 for risk aversions log 2 7 in the range [—5,5] and 
maturities T in the range [0.1, 1] years. Figure 4 shows the dollar amount of the hedge in 
excess of the Merton strategy for the volatility put as a function of Y and T for 7 = 1. 
These experiments showed in particular that there is no significant dependence of these 
values on 7 in the range [0, 10]. We can also observe how, for longer maturities, the mean 
reverting character of Y t leads prices which are practically independent of Y . 

The present model has the somewhat surprising property that unbounded claims such 
as forward contracts and call options have an expected utility of negative infinity. For 
that reason, we restrict ourselves to claims for which the payoff is bounded from above, 
such as the difference of calls. Figure 5 shows the indifference price for a difference of 
volatility calls {a\ - K ± ) + - (erf, - K 2 ) + with K x = 0.15, K 2 = 0.3 for maturities T in the 
range [0.1, 1] years and risk aversion 7 = 1. 

In figure 6 we select one sample path for the squared volatility process Y t = of and the 
indifference price for the volatility put with K = 0.15, T = 1,7 = 1. We then show the 
difference between the amount invested in the stock according to the hedging portfolio and 
the amount invested in the stock following Merton's strategy. This plot shows clearly that 
long before maturity, the hedger can ignore the contingent claim, but that they become 
increasingly anxious closer to maturity. Further experiments show again that these values 
are rather insensitive to the choice of 7 in the range [0, 10]. 
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8 Conclusions 



The main goal of this paper has been to find a financially viable class of stochastic 
volatility models and related derivative securities which lead to tractable forms for the 
utility based indifference price and the corresponding optimal hedging strategies. As is 
well known, the indifference price is wealth independent only for exponential utility, so 
we have adopted this as the preferred utility. The indifference price and its linearization 
the Davis price are motivated by fundamental economic principles applicable to general 
incomplete markets, which makes these questions of theoretical interest and of use to 
practitioners using models with either stochastic volatility or jumps. 

We have found that "close to closed form solutions" can be derived very elegantly for 
contingent claims which depend solely on spot volatility at maturity. A class of stochastic 
volatility models for which this is possible was found which is as rich as but inequivalent 
to the original Heston model. Contingent claims on spot volatility have been studied 
before in the literature, but always with an implied assumption about the market price of 
risk. Our treatment is the first to use the market price of risk derived from first principles 
via the utility function. A number of the interesting properties of volatility claims have 
been explored by other authors (such as in Griinbichler and Longstaff 1996), and these 
effects can be easily observed in our formulation. 

Our formulas can be efficiently computed using the fast Fourier transform: the result- 
ing algorithm yields hundreds of prices per second on a desktop PC. This speed gives us 
the capability to explore the indifference price, optimal hedging and sensitivities and the 
market price of risk for a complete range of financially relevant model parameters includ- 
ing risk aversion, mean rates of return, mean reversion rates, correlation parameter, mean 
volatility and volatility of volatility. We have shown a number of the possibilities in this 
paper, but many detailed properties remain to be explored in further papers. 
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Amongst the important questions still to be addressed we mention two. The first 
is to quantify the performance of our optimally hedged portfolios relative to unhedged 
portfolios for a variety of volatility claims. Since the volatility claims we analyze are 
traditionally considered to be "unhedgable" , it is of interest to use quantities such as 
value-at-risk to measure portfolio performance. We expect to observe performance which 
is directly related to the correlation between volatility and stock price. A second question 
is to compare the exponential utility based market price of risk we compute with the 
forms assumed in the literature by other researchers. We hope to show that our computed 
market price of risk gives reasonable consistency with observed option price data. 
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